function u = fn_u(c,gamma)

u = (c.^(1-gamma))./(1-gamma);

end